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SECTION 1 


INTRODUCTION 


The construction and maintenance of the primary reflectors of future large 
space telescopes will demand major advances in the technology of space struc- 
tures and materials. For example, the so-called Large Deployable Reflector 
(LDR) , described in Reference 1, Is planned to have an aperture of 20 meters 
and a wavefront accuracy of about five micrometers. The resulting ratio of 
allowable structural error to aperture of 10“ 7 is much smaller than present 
antenna systems and is surpassed only by non-deployable advanced optical 
instruments such as the Hubble Space Telescope. 

Large reflectors such as the 20-meter LDR cannot, of course, be transported to 
space in their final configuration. In addition, because of the requirement 
that the surface of the reflector be mirror-like and very stable, the con- 
struction by purely deployable techniques is difficult. One such approach, 
described in Reference 2, and illustrated in Figure 1, utilizes a deployment 
canister containing interconnected modules of reflector panels and attached 
stowed truss segments. The canister walks Itself around the perimeter of the 
reflector, deploying the modules sequentially and attaching them to the grow- 
ing structure. Such an intelligent deployment mechanism would require exten- 
sive development. 

The alternative to pure deployment is to use assembly in space. A possible 
approach would be to attach premanufactured reflector panels to a supporting 
stiff truss structure. While the truss structure could also be erected, the 
total on-orbit assembly work would be greatly reduced by deploying it. The 
recent invention of the Pactruss concept, described in Reference 3, with its 
strongly self-synchronized deployment, makes this combined approach particu- 
larly attractive. 

The application of the Pactruss deployable structure to the Space Station pri- 
mary structure is described in References 4 and 5. The concept is seen to 
provide a realistic alternative to assembly in enabling the construction of 
large trusses from small stowed volumes. In a similar way, the usefulness of 
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the concept to precision segmented reflectors can be shown. For high 
concentration-ratio solar collectors, for example, the Pactruss concept is 
shown in Reference 6 to provide an efficient precise backup to the reflector 
panels. 

Because of its apparent ease of synchronous deployment, as well as its excel- 
lent deployed dimensional stability, the Pactruss concept has been identified 
as the primary candidate for the deployable truss for LDR (see Reference 1). 
Some primary issues, however, need to be Investigated. One is whether the 
Pactruss design with its stringent geometric requirements can indeed be 
applied to a doubly-curved surface such as the parabolic reflector of LDR. 
Another is how to modify the design so that it can be stowed around, and 
attached to, a central part of the spacecraft for launch. In addition, quan- 
titative evaluations of the weight and stiffness of the resulting deployed 
truss structure are needed. The investigations are described herein. 
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SECTION 2 


GEOMETRY OF HYBRID PACTRUSS 


Pactruss deployable structures stow In a very compact form. As shown schemat- 
ically in Figure 2 and as a conceptual model in Figure 3, the structure is 
composed of verticals connected by surface members (longerons) and core mem- 
bers (diagonals). The verticals move up or down during stowage, the sense 
alternating in two of the three triangular directions over the planform. 
Alternating verticals are connected by non-folding longerons. In the third 
planform direction, adjacent verticals move in the same sense. The longerons 
in this direction must fold in order to allow the ends to move together. The 
verticals and the non-folding longerons form sets of parallelograms that 
become thin when stowed. In the stowed configuration, all members are 
vertical, with the foiding longerons and diagonals occupying spaces provided 
by finite hinge offsets. The stowage is one of the type called "double fold" 
because the planform is packaged in both directions. 

For a configuration such as the LDR, Figure 4, in which the deployed truss 
surrounds a center body, double-fold stowage presents a problem— the central 
opening occupied by the center body also packages during stowage. Thus, the 
stowed truss must be packaged separately from the center body and mounted to 
it after deployment. A better approach would be to have the truss stow around 
the center body. It could then be attached and any utility connections could 
be made and tested prior to launch. 

The desirability of stowage around the center body motivated the creation of 
the hybrid Pactruss concept shown in F'igure 5. The hybrid is made up of 
single-fold beams separating the planform Into areas which are filled with 
double-fold Pactruss. The directions of the non-folding Pactruss longerons 
are taken to be parallel to the beams on the radial boundaries. On stowage, 
the hybrid truss shrinks to a "sleeve" around the center body as shown in 
Figure 6. The single-fold beams accordion-fold to flat packs and the Pactruss 
segments occupy the small triangular regions at the corners. 

The width of the single-fold beams should be selected so that the central 
opening fits properly around the center body. If possible, the width also 
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should be selected to provide good locations for the attachment of reflector 
panels. A possible arrangement of panels is shown by the heavy lines in 
Figure 7. Here three corners of each panel are attached to the truss node 
lying under them. In the single-fold regions, the attachments are located 
differently than in the double-fold regions. 

The truss configuration in Figures 5 through 7 is divided into six congruent 
sections. It is entitled "Sixpac" and is shown in perspective in Figure 8. 
Another possible arrangement would be the "Tripac" shown in Figure 9. Note 
that the beams are considerably wider in this particular version in order to 
interface with larger panels attached as shown by the heavy dots in the 
figure. Note that each panel has its own three attachment points instead of 
having to share them with adjacent panels. 

Of course, other arrangements are possible. The two shown serve to demon- 
strate that pre-attachment of the truss to a center body is feasible. 
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SECTION 3 


GEOMETRIC REQUIREMENTS 
OF THE PACTRU3S CONCEPT 


The basic geometrical requirements for all trusses with slender members is 
that their centerlines must pass through a common paint at the node to which 
they are attached. 

The primary geometrical requirement pertaining to the Pactruss concept is that 
the distance between connected nodal pairs be maintained throughout deploy- 
ment. Otherwise, the misfit would generate large axial loads in the members. 
For the folding longerons and diagonals, this requirement can be easily met by 
appropriately locating the knee hinge in each member. The non-fold mg 
members, on the other hand, must span the correct distance exactly. 

The secondary geometrical requirement is that the hinges between members be 
located and oriented so that the amount of bending and torsion in each strut 
is small. Optimally, the entire deployment should be strafim-free, except for 
those strains Induced by hinge moments (either driving or frictional) or 
external loads.' 

Both of these requirements are complicated by the fact that the struts, 
although slender, have finite cross-sectional dimensions and require space in 
the package. Therefore, the hinges between the members must be offset from 
the center of the nodal points of the truss by as much as three member 
diameters. Preferably, the hinges at member ends should be located on their 
centerline so that no upsetting moments are generated by axial forces in the 
members of the deployed truss. 

When the Pactruss forms a flat surface, with regularly spaced modules, the 
design of the Pactruss joints is straightforward. All the non-folding 
longeTdrrs—a-p* made the same length, with the same total hinge offset at the 
two ends, producing a deployed structure with uniform distances between the 
node points along the directions of the non-folding longerons. The hinge 
lines are oriented in the horizontal plane, perpendicular to the member 
centerline. During deployment, the projection' of each node in the surface 
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plane moves directly and proportionally away from the projection of its 
neighbors. The deployment proceeds smoothly, with the non-folding longerons 
at a uniformly decreasing slope with respect to the horizontal plane. 

For a doubly curved surface, on the other hand, the proper design is not so 
immediately evident. Careful examination, however, shows that at least one 
solution exists that satisfies the primary requirement of correct distances 
between nodes. Furthermore, the solution also yields hinge locations and ori- 
entations that satisfy the secondary requirement for the conditions of full 
deployment and full stowage. Whether the solution also allows strain-free 
deployment is by no means evident; the evaluation of deployment' straining 
requires detailed analysis. 

The design solution is an extension of that for the flat truss. First, the 
distances between adjacent deployed nodes along non-folding longerons must be 
a constant value. Second, the sum of the hinge offsets at the ends of each 
non-folding longeron must be the same. Third, the hinge lines must lie in the 
horizontal plane and be perpendicular to the longeron centerline. Note that 
these are the same as the rules used for the flat truss, except that the off- 
set lengths can be different for the two directions in the flat truss— not so 
with the doubly curved truss. 

That the solution yields a deployed and stowed strain-free Pactruss can be 
seen from the following reasoning: 

Let the distance between nodes connected by non-folding longerons be l and the 
total of the hinge offsets at the two ends be e. Let the upper-surface points 
of the deployed truss be located at (x n> y n , z n ) n = 1,2,3...,N with 

(X, - Xj) 2 ♦ (y, - y/ = (2, - Zj ) 2 - i 2 

for all i , j pairs connected by non-folding longerons. Assume (without loss of 
gerrcraiityy that the origin of the coordinate system is attached to one cf the 
nodes which does not move upward during stowage. Let the fixed node number be 
m. Then the coordinates of the n^ h node when stowed can be seen to be 




where z u * 2 - e, for upwards stowing nodes 
* 0, otherwise 

One way to visualize this stowage Is to assume that the non-folding longerons 
are removed and the hinge lines on either -end are merged. The resulting nodal 
surface is then the same as the deployed one, except that it is shrunken by 
the factor e/2. Inserting the longerons then raises the upward stowing nodes 
by the distance 2 - e. 

Examination of the hinge lines shows that they have the same orientation 
whether deployed or stowed; the verticals have no rotation. The longerons 
thus do not bend or twist. 

In order that the nodes lie on the desired paraboloidal surface and obey the 
internodal distance requirement, a computational algorithm must be developed 
as follows: 

Let (x^ y^, Zj) and (x 2> y 2> z 2 ) be previously located points. The location 
of the next point (x, y, z) must satisfy the following equations: 

(x - x 2 ) 2 + (y - y x ) 2 + (z - z^ 2 = 2 2 ' 


(x - x 2 ) 2 + (y - y 2 ) 2 + (z - z 2 ) 2 = 2 2 



7 



where F is the focal length of the parabaloid. 

The easiest way to solve for x, y, and z is to assume a value for z, solve the 
first two equations for x and y, and determine a new value for z. Repeating 
the process a finite number of times (say, 50) guarantees a very accurate 
answer. 

Before determining the Pactruss nodes, the nodes for the beam lying along the 
x axis are found by using only one of the distance equations in addition to 
the paraboloid equation, with the assumption that y is the half-width of the 
beam. Thus for the beam of width b, 


(x - x.) 4 + (z - z.r = V 


b 2 x 2 
Z 55 Z o + I6F + 4F 


This is also solved most easily by iteration. 

After the points on the upper surface are obtained, those for the lower sur- 
face are determined by subtracting the constant depth h from the z coordi- 
nates. 

In order to determine the points, computer programs were written in the C 
language. The source code is included in the Appendix for GEN6PAC.C and 
GEN3PAC.C for generating the Sixpac and Tripac geometries, respectively. 
Inputs to these programs are: 

N, the number of rings 
4, the longeron length 
h, the truss depth 
b, the beam width 
F, the focal length 

Output is a file, TRUSSNAME.DAT, listing the nodal coordinates. 

Also included in the Appendix are source files GEN6STR.C and GEN3STR.C, which 
read the file TRUSSNAME.DAT and generate a file, TRUSSNAME.DTA, which contains 
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the nodal coordinates followed by a listing of node endpoints and strut type 
for each of the members. 

An example of an output file for a two-ring Slxpac truss is included in 
Table 1. Note that there are 72 nodes (36 on each surface), 258 members, and 
12 types of members. An additional point is included at the focal point. For 
this case, i = 2.0, h = 2.0, b = F = 5.0. 
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SECTION 4 


STRUCTURAL PERFORMANCE 


In order to obtain quantitative estimates of the ability of the Pactruss 
structure to provide precise positioning and support to reflector panels, the 
five-ring Sixpac truss was analyzed using the MSC/pal 2* finite element 
analysis. The truss is shown in Figure 8 and has the following dimensions: 


F = 10.0 m 

A = 2.0 m 

h = 2.0 m 

b * 1.732 m 


The necessary input file to the MSC/pal 2 program was generated by the com- 
puter program SEE2PAL.C .included in the appendix, which reads the 

TRUSSNAME.DTA file previously prepared, and obtains structural property infor- 
mation from the data file TRUSSNAME.PRP. This latter file lists the strut 

area, modulus, and so forth, for each strut type. 


4.1 STATIC PERFORMANCE 

Static structural analysis was performed for two cases: One-g loading In the 

negative z-direction, and an angular acceleration of 0.001 radians/sec 2 around 

the y-axis. The first case was selected to determine the expected deflections 

during ground testing. The second case represents a severe case of loading 

due to steering control while in operational orbit. (Earlier estimates given 

in Reference 2 reveal that the maximum operational angular acceleration is of 
-4 2 

the order of 10 radians/sec .) In both cases, the truss is assumed to be 
mounted at the central opening. 

Results of the static analyses, given in Figure 10, pertain to both the bare 
truss alone and to the truss loaded with reflector panels weighing 10 kg/m 2 . 


*MSC/pal 2 is a trademark of the MacNeal-Schwendler Corporation. 
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The mass of the structure itself included a joint weight equal to the weight 
of the tubular struts. 

Examination of the deformation shapes showed that most of the deformation 
comes from straining of the inner rings. Accordingly, a tapered-stiffness 
design was developed which was much lighter than a uniform stiffness truss of 
equal overall stiffness. The tapered design was chosen for further study. 

The tapered design deflects a maximum of 200 micrometers at the rim due to its 
own weight. The accuracy requirement that has been placed on the truss is 
that its attachment nodes should be precise to 100 micrometers without on- 
orbit adjustment. The fact that the gravity sag is only twice this much gives 
confidence that a combination of ground testing and analysis will be able to 
guarantee the required accuracy when the truss is deployed on-orbit. 

The structural weight, including joints, of the tapered design is 819 kg. The 
truss supports panels weighing 3730 kg and is strong enough to carry them in a 

one-g field. Of more import is the fact that the maximum deflection at the 

2 

rim caused by the angular acceleration of 1.0 milliradian/sec is only 3.5 
micrometers. 

The distortion due to operational loading Is better described in terms of the 
rms deviation from a best-fit paraboloidal surface. The source code of the 
program BESTFIT.C is included in the Appendix. This program reads the files 
output by MSC/pal 2 and determines the deviation from and the pointing error 
of the best-fit paraboloid. For the t'wo-meter deep truss, the rms deviation 
is only 0.17 micrometer; the pointing error is 2 x 10" 7 radians or 0.04 arc 
sec. These errors are well within the budgets of ~3 micrometers and 0.1 arc 
sec established for the LOR in Reference 1. The fact that the Pactruss is 
stiff enough to prevent unacceptable deflections due to operational loads 
means that the control system needed to correct the optical path can be of the 
low band pass type and interaction between the structure and the control sys- 
tem can be readily handled. 

4.2 EFFECTS OF TRUSS DEPTH 

Some influences of truss depth on structural performance are shown in Figure 
11. The distortions increase rapidly for truss depths less than two meters. 
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Larger depths reduce the distortion at the cost of increased structural mass. 
Of more concern is the fact that the longer struts have a lower lateral vibra- 
tion frequency, which is about 22 Hz for a two-meter depth. Earlier studies 
of sensitivity of surface accuracy to manufacturing imperfections (Reference 
6) indicate that depths between one and three longeron lengths are best. All 
these results point to the desirability of a truss depth equal to or moder- 
ately greater than a longeron length. 

4.3 VIBRATION PERFORMANCE 

As an additional indication of the stiffness of the Pactruss, the dynamic 
analysis capability of MSC/pal 2 was used to obtain the natural vibration 
modes and frequencies for the tapered design. The structure includes the 
truss and the panels and is mounted at the 12 node points at the central 
opening. The results are shown in Figure 12 for the first 20 modes. The fun- 
damental frequency is over 10 Hz, a high figure indeed for such a large 
structure. The modal density is high, with the first five modes having fre- 
quencies within 15 percent of each other. 
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SECTION 5 


DEPLOYMENT ANALYSIS 


In order to obtain some insight Into the behavior of a Pactruss structure with 
a doubly curved surface, the sector shown darkened in Figure 13 was analyzed 
with an Astro Aerospace Corporation proprietary deployment program called 
ASTRAN. This computer analysis program was written during 1987 and has the 
following characteristics: 

• High fidelity, efficient analysis of deformations and loads 
encountered by flexible truss structures during deployment. 

• Large displacements and rotations, small distortions, quasi- 
static. 

• Structural elements are clusters (hinge bodies) and struts 
connected by hinges. 

• Can apply external loads, hinge moments, constraints, and pre- 
scribed hinge angles. 

• Degrees of freedom are six for each cluster and one for each 
hinge. 

• Written in C language; presently running on PC-family com- 
puters. 

• ASTRAN is interactive. The analysis can be Interrupted in 
order to examine aspects of the behavior and to change updat- 
ing procedures in order to help convergence. 

. • The user develops a good "feel" of how the structure behaves. 

The similarity to conducting a test is uncanny. 


The method has been applied to a two-bay segment of articulated Astromast con- 
sisting of 15 clusters, 33 hinges and 27 struts. The results obtained com- 
pared well with experimental results. It was also used to predict the 
stability of deployment of Z-Beam, which was being developed for the now- 
defunct COFS Mast Flight Experiment program. Results showed successful stow- 
age even with a large imbalance in hinge friction moment. 

The effect of the surrounding structure on the segment shown In Figure 13 was 
simulated by constraining the nodes along the beams to have a zero component 
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of displacement normal to the plane of the beams. The deployment is assumed 
to be driven at the two hinges joining the lower non-folding longerons to the 
inner cluster. The outer bays were allowed to deploy freely. 

The analysis revealed a serious deployment defect in the Sixpac design. The 
inner bay deployed easily. The outer bays also deployed well for the initial 
stages, exhibiting no member straining, but as the inner bays reached full 
deployment, the outer bays "hung up" in the partially deployed condition shown 
In Figure 13. The difficulty arises from the fact that the folding longerons 
straighten prematurely and lock the outer bays, whose longerons should swing 
through the horizontal from the stowed downward sloping condition to the 
deployable upward sloping condition. 

Numerous trials were made attempting to get the outer bays to pop through 
before the inner bay was fully deployed. Only by providing large asymmetri- 
cal drive moments to the outer bays could full deployment be achieved. 

A model based on the Tripac design was similarly analyzed. Here the results 
were more promising. It is possible to achieve full deployment by driving the 
bays whose longerons must pass through the horizontal. Successful deployment, 
however, will probably require driving at many locations, thereby destroying 
the main advantage of the Pactruss concept. 

One possible, remedy is to modify the Pactruss concept somewhat in a manner 
shown in Figure 14. In this approach, the situation in which a bay's longe- 
rons must pass through the horizontal is avoided by requiring such bays to be 
horizontal when deployed. Extensions, with braces if necessary, provide the 
required mounting point at the paraboloidal surface. Bays for which the 
longerons move through less than 90 degrees are designed to be steeper, 
thereby returning to the paraboloidal surface. 

The results for the deployment analysis of a sector of the modified truss 
shown darkened in Figure 14 and called XGTK are given In Figure 15. The 
deployment angle of the outer bays Is shown as a function of the drive angle 
at the inner bay. The results for XGTK show smooth progress to full deploy- 
ment even with no outer helping hinge moment, M h . In comparison, the results 
for the unmodified design, identified as XGTEST, show smooth deployment until 
nearly full deployment, where the structure jumps to the undesirable result. 
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Clearly, fuller models of the deploying truss must be investigated in evaluat- 
ing the deployment behavior. The effects of off nominal drive and hinge con- 
ditions and the effects of external forces need to be examined. These 
preliminary results indicate that the modified Pactruss with extensions will 
deploy successfully. 
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SECTION 6 


CONCLUDING REMARKS 


This investigation has shown that it is possible to design a hybrid Pactruss 
structure capable of supporting reflector panels to form a precise, doubly 
curved (paraboloidal) reflector surface. The truss can be stowed around a 
central body to which it is attached for launch and orbital deployment. The 
deployed truss is very stiff, stiff enough to resist operational loads with 
sub-micrometer rms distortions and to limit gravity sag to amounts small 
enough to enable ground testing and measurement of the surface. 

Preliminary deployment analyses indicated that the Pactruss, modified with 
extensions to avoid lockup, can be deployed in a strain-free manner. 
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Figure 1 . Sequentially deployable precision reflector. 
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STOWED 


Figure 2. LDR Paclruss deployment scheme. 
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BLACK AND WHITE PHOTOGRAPH 



(a) Cutaway perspective 



(b) Cross section 


Figure 4. Application of Pactruss to the support structure 
for the primary mirror of the JPL low-cost LDR. 
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Figure 5. Hybrid Pactruss concept. 
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Figure 6. Stowed and deployed hybrid Pactruss. 



Figure 7. Hexagonal panels attached to hybrid Pactruss. 



Figure 8. Sixpac support truss. 



Figure 9. Tripac hybrid Pactruss support structure with panels. 
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CONFIGURATION STRUCTURAL PANEL MASS MAXIMUM DISPLACEMENT (nm) 

5-Ring, 2-Meter Sixpac MASS (kg) (10 kg/m 2 ) (kg) (1 g) lx 10~ 3 rad/sec 2 ) 
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Figure 10. Static performance for the hybrid Pactruss (E = 30.0 x 10 10 N/m 2 , Kjoint = 2). 
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Figure 1 1 . Effect of truss depth. 


MODE 


MODE SHAPE 



1 

1.03718E+01 

N - 1, ROCKING A80UT XX 

2 

1.03718E+01 

N * 1, ROCKING ABOUT YY 

3 

1. 17236E+01 

N * 0, ROTATION ABOUT ZZ 

4 

1.19297E+01 

N » 2, NODE LINES AT -45,45 DEG 

5 

1.19297E+01 

N - 2, NODE LINES AT 15,105 DEG 

6 

1.37824E+01 

N « 0, AXIAL PUMPING 

7 

1.59909E+01 

N * 3, NODE LINES AT -30,30,90 DEG 

8 

1.95641E+01 

N - 3, NODE LINES AT 0,60,120 DEG 

9 

2.47072E+01 

N . 4, NODE LINES AT 0,45,90,135 DEG 

10 

2.47072E+01 

N - 4, MAXIMA AT 0,45,90,135 DEG 

11 

2.84482E+01 

N - 5, MAXIMA AT -54,18,90,162,234 DEG 

12 

2.84482E+01 

N * 5, NODE LINES AT -54,18, ETC. DEG 

13 

3.26987E+01 

N - 5, MAXIMA AT -54,18, ETC. DEG 

14 

3.26987E+01 

N - 5, NODE LINES AT -54,18, ETC. DEG 

15 

3.52508E+01 

N - 6, MAXIMA AT 0,60, ETC. DEG 

16 

3.91162E+01 

N * 6, NODE LINES AT 0,60, ETC. DEG 

17 

4.06708E+01 

N - 3, SECOND RADIAL MODE; NODES AT 0,60,... 

18 

4.06708E+01 

N * 3, ?? 

19 

4.21089E+01 

N * 4, ?? 

20 

4.21089E+01 

N * 2, SECOND RADIAL MODE; MAX AT 0,90,... 


Figure 12. Vibration frequencies for center-mounted Slxpac truss 
with tapered truss stiffness and 10 kg/m 2 panels. 
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Figure 13. XGTEST, a model for investigating deployment 
of the sixpac configuration. 
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Figure 15. Results of outer-bay deployment for two models. 



APPENDIX A 

Example Files 

GEN6PAC.C 
* GEN3PAC.C 
GEN6STR.C 
GEN3STR.C 
TWORING.DTA 
SEE2PAL.C 
BESTFIT.C 


A-l 


*• * 


I 


/* GEN6PAC.C - Program for generating the basic geometry of a hybrid Pactruss 

* with six radial single-fold beams bounding the six sectors of 

* double-fold Pactruss. The truss nodes are required to lie on 

* the paraboloidal surface 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

.* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

♦ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

*/ 


z = r*r/(4*F) 

where F is the focal length and r is the radius 


r = sqrt(x*x + y*y) 

The truss is composed of M rings. The number of nodes in each 
sector is (M + 1)*(M + 2)/2. The length of the non-folding 
surface struts is 1, and the width of the single-fold beams is b. 
The truss height is h. 

Although this analysis does not explicitly consider the effects 
of joint offsets, the geometry is exact for the case in which the 
length 1 is divided between offsets, totalling 3.5 tube diameters 
times the secant of the surface slope at the rim, and the strut 
itself. The stowed geometry then can be obtained by uniformly 
shrinking the planform and appropriately adjusting the z-coordi- 
nate of each node. Note that the deployment mechanics may 
necessarily involve some straining of the members as a result of 
the double curvature of the nodal surface, not as a result of 
finite offsets. 

The nodes are numbered consecutively within each sector, starting 
at the first sector with 1, 2, 4, 7, 11, 16, etc. outward along 
the beam boundary with intermediate numbering counterclockwise 
along each ring. Nodes in succeeding sectors are numbered in the 
same order consecutively. The entire upper surface is numbered 
first, then the lower, and finally any accessory nodes are as- 
signed numbers. 

The focal point is assigned the number 0. 

The results are stored in ASCII form in the file XYZ.DAT where 
XYZ is entered in the command tail. The default is SIXPAC. 

The first entry is the parameter values. The second line contains 
the total number of nodes in the truss. The following lines are 
the coordinates of each node in x,y,z form listed in order. 

The input parameters M, F, 1, b, and h are obtained from* the file 
XYZ.PHM, which is a one-line ASCII file with the values, in order, 
separated by commas. The default values for SIXPAC are 6, 10.0, 
2.0, 1.154701, and 2.0. 

John M. Hedgepeth 8/26/87 

Altered to allow F = infinity by inputting F < 0 9/16/87 

Updated to MS C 5.0 12/29/87 


♦include <scitech.h> 


int M = 5; 
double F = 10. ; 
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double 1 = 2. ; 

double b = 1.154701; 
double h = 2.; 

void initbeam( ) ,newpoint( ) , fillsec( ) ; 

main(argc,argv) 
int argc; 
char *argv [ ] ; 

{ 

int i,j,n,chr,N; 

static double x,y,z,coord[100] [3] ,cosv,sinv; 
static char inname [ 50 ] ,outname( 50] ; 

FILE *infile,*outfile; 

strcpy(outname, "SIXPAC.DAT") ; 

if(argc > 1) { 

strcpy(outname,argv[l] ) ; 
strcat(outname, ".DAT") ; 
strcpy( inname, argv[l] ) ; 
strcat ( inname , " ; PRM" ) ; 
strupr( inname) ; 

if ( ( infile = fopen(inname,"r")) == NULL) { 

printf("\nCannot open %s. Aborted. \n", inname) ; 
exit(l) ; 

} 

if (fscanf (infile," Xd Xlf *lf *lf *lf", &M, &F, &1, &b, &h) != 5) { 
printf("\nCannot read input parameters. Aborted. \n" ) ; 
exit(l); 

} 

fclose( infile) ; 

if(F == 0) 

F = -10000.5 

> 

N = ((M + 1)*(M + 2) )/2; 

if((outfile = fopen(outname, "r") ) != NULL) { 
fflush(stdin) ; 

printf("\nFile Xa exists. Do you want to write over it? (Y/N) <N>", 

outname) 

if(toupper(getch()) ! = * Y* ) 
exit(l) ; 

fclose(outfile) ; 

> 

outfile = fopen(outname,"wt"); 

fprintf (outfile, " ^d, %lf, Xlf, )tlf, %lf\n",M, F, 1, b, h) ; 
fprintf (outfile," %d\n",12*N +1); 
x = 0; 

fprintf (outfile, " %.81e, *.81e, %.81e\n",x,x,F) ; 


initbeam( coord) ; 
fillsec(coord) ; 


cosv = cos( (double) i*PI/3) ; 
sinv = sin( (double) i*PI/3) ; 

for(n=l; n<= N; n++) { 

x = coord[n] [0]*cosv - coord[n] [l]*sinv; 
y = coordfn] [oj*sinv + coordfn] fl]*cosv; 
z = coordfn] [2] ; 

fprintf (outfile, " *.81e, %.81e, %.81e\n”,x,y,z) 

} 

} 


for(i=0; i<6; i++) { 

cosv = cos( (double) i*PI/3) ; 
sinv = sin ( (double) i*PI/3) ; 

for(n=l; n<= N; n++) { 

x = coord [n] [0]*cosv - coord[n] [l]*sinv; 
y = coordfn] fO]*sinv + coordfn] fl]*cosv; 
z = coordfn] f2] - h; 

fprintf (outfile, " X.81e, %.81e, *.81e\n",x,y,z) 

> 

} 

fclose(outfile) ; 

printf ("\nFile *s successfully written. \n", outname) ; 


void initbeam( coord) 
double coord[100] [3] ; 

{ 

int i,n, j; 

double x,y,z, xold, ksi,sqrt3, eps; 

sqrt3 =3.; 

sqrt3 = sqrt(sqrt3); 

x = sqrt3*b/2. ; 
y = b/2. ; 

coord[l][0] = x; 
coordfljfl] = y; 
coord[l] [2] = b*b/(4.*F); 
if (F < 0) 

coordfl] [2] = 0; 

for(i=l; i<=M; i++) { 

n = (i*(i + l))/2 + 1; 
if (F < 0) { 
x += 1; 

coord[n] [2] = 0; 

} 

else { 

' xold = x; 

eps = 1/2. ; 

for(j=0; j<50; j++,eps/=2.) { 
ksi = x + eps - xold; 
z = ksi*(ksi + 2.*xold)/(4.*F) ; 
if(ksi*ksi + z*z <= 1*1) 
x += eps; 



coord[n][2] = (x*x + y*y)/(4.*F); 

} 

coord[n] [0] = x; 
coord(n][l] = y; 

} 

} 


void fillsec( coord) 
double coord [100] [3] ; 

{ 

int i,n,p,q,r; 

double x,y,z, temp, cos v,sinv; 

temp = PI/3. ; 
cosv = cos (temp); 
sinv = sin(temp); 

/* First fill in the other radial boundary 
*/ 

for(i=l; i<=M; i++) { 

n = ((i + l)*(i + 2))/2; 
x = coord[n - i] [0] - b*sinv; 

coordfn] [0] = x*cosv + b*sinv; 
coord[n][l] = x*sinv + b*cosv; 
coord[n] [2] = coord[n - i][2]; 

} 

/* Then fill horizontal rows 
*/ 

for(i=l; i<=M-l; i++) 

for(n=i+l; n<=M; n++) { 
p = ((n - l)*n)/2 + i; 

q = p +. l; 

r = (n*(n + l))/2 + i + 1; 
newpoint(p,q,r, coord) ; 

} 

} 


void newpoint(p,q,r, coord) 
int p,q,r; 

double coord[100] [3] ; 

{ 

int i,n; 

double x,y,z,zold,Id}ar[3] ,delR[3] ,rhosq,d, temp.eps, mult; 

for(i=0,rhosq=0; i<3; i++) { 

Hbar[i] = (coord[p][i] + coord[q] [i] )/2; 
delRfi] = coord[p][i] - coord(q][i]; 
rhosq += delR[i]*delR[i] ; 

} 

d = rhosq - delR[2]*delR[2] ; 

zold = (Rbar[0]*Rbar[0] + Rbar[l]*Rbar[l] )/(4.*F) ; 
eps = zold/16; 

for(i=0,mult=l. ; i<50; i++,eps*=mult) { 
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z = zold + eps - Rbar[2]; 
if (F < 0) { 

z = zold = eps = 0; 
i = 49; 

} 

temp = 1*1 - rhosq*(0.25 + z*z/d) ; 

temp /= d; 

temp = sqrt(temp); 

x = Rbar[0] - z*delR[0]*delR[2]/d - delR[l]*temp; 
y = Rbar[l] - z*delR[l]*delR[2]/d + delR[0]*temp; 

if(x*x + y*y >= 4.*F*(zold + eps)) 
zold += eps; 
else mult = 0.5; 


} 

coord[r][0] = x; 
coord[r][lj = y; 
coord[r][2] = zold; 
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/* GEN3PAC.C - Program for generating the basic geometry of a hybrid Pactruss 
with thr.ee radial single-fold beams bounding the three sectors of 
double- fold Pactruss. The truss nodes are required to lie on 
the paraboloidal surface 

z = r*r/(4*F) 

where F is the focal length and r is the radius 


r = sqrt(x*x + y*y) 

The truss is composed of M rings. The number of nodes in each 
sector is M*(M +3) +1. The length of the non-folding 
surface struts is 1, and the width of the single-fold beams is b. 
The truss height is h. 

Although this analysis does not explicitly consider the effects 
of joint offsets, the geometry is exact for the case in which the 
length 1 is divided between offsets, totalling 3.5 tube diameters 
times the secant of the surface slope at the rim, and the strut 
itself. The stowed geometry then can be obtained by uniformly 
shrinking the planform and appropriately adjusting the z-coordi- 
nate of each node. Note that the deployment mechanics will 
necessarily involve some straining of the members as a result of 
the double curvature of the nodal surface, not as a result of 
finite offsets. 

The nodes are numbered consecutively within each sector, starting 
at the first sector with 1, 3, 7, 13, 21, etc. outward along 
the boundary of the beam running out along the x axis . Along the 
other beam the points are numbered 1, 4, 7, 14, 22, etc. The 
other points in the first sextant is numbered with odd numbers 
increasing toward the center. The even numbers are assigned to 
the symmetrically located points in the second sextant. Points are 
assigned to the midpanel points of the beam for possible use; they 
are even numbered 2, 6, 12, 20, 30, etc. Nodes in succeeding 
sectors are numbered in the same order consecutively. The entire 
upper surface is numbered first, then the lower, and finally any 
accessory nodes are assigned numbers. 

The focal point is assigned the number 0. 

The results are stored in ASCII form in the. file XYZ.DAT where 
XYZ is entered in the command tail. The default is SIXPACK. 

The first entry is the parameter values. The second line contains 
the total number of nodes in the truss. The following lines are 
the coordinates of each node in x,y,z form listed in order. 

The input parameters M, F, 1, b, and h are obtained from the file - 
XYZ.PHM. The default values for TRIPAC are 5, 10.0, 2.0, 3.4641, 
and 2.0. 

John M. Hedgepeth 9/6/87 

Updated to MS C 5.0 12/29/87 


♦include <scitech.h> 
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int M = 5; 
double F = 10.; 
double 1 = 2. ; 

double b = 3.4641; 
double h = 2. ; 

void initbeam() ,sympoint() ,newpoint( ) , fillsec( ) ; 

main(argc,argv) 
int argc; 
char *argv [ ] ; 

{ 

int i,j,n,chr,N; 

double x,y,z,coord[100] [3] ,cosv,sinv; 
char inname[50] ,outname[50] ; 

FILE *infile,*outfile; 

strcpy(outname, "TRIPAC.DAT") ; 

if(argc > 1) { 

strcpy (outname , argv [ 1 ] ) ; 
strcat ( out name , ” . DAT” ) ; 
strcpy( inname, argv[l] ) ; 
s t rcat ( inname , " . PRM" ) ; 
strupr(inname) ; 

if ( ( infile = f open ( inname, ”r”)) == NULL) { 

printf(”\nCannot open %a. Aborted. \n”, inname) ; 
exit(l) ; 

} 

if (fscanf( infile," *d *lf %lf %lf *lf”, &M, & F, &1, &b, &h) != 5) { 
printf("\nCannot read input parameters. Aborted. \n") ; 
exit(l) ; 

} 

f close ( infile) ; 

} 

N = M*(M + 3) + 1; 

if((outfile = fopen(outname,"r”)) != NULL) { 
fflush(stdin) ; 

printf("\nFile %s exists. Do you want to write over it? (Y/N) <N>", 

outname) ; 

if(toupper(getch()) != *Y*) 
exit(l); 

fclose(outfile) ; 

} 

outfile = fopen( outname, "wt") ; 

fprintf (outfile, ” *d, %lf, Xlf, S;lf, S;lf\n",M, F, 1, b, h); 
fprintf (outfile," *d\n",6*N + 1); 
x = 0; 

fprintf (outfile,” *.81e, *.81e, S;.81e\n”,x,x,F) ; 


initbeam(coord) ; 
fillsec( coord) ; 
for(i=0; i<3; i++) { 
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cosv = cos ( (double) i*PI*2./3. ) ; 
sinv = sin( (double) i*PI*2./3. ) ; 

for(n=l; n<= N; n++) { 

x = coordfn] [0]*cosv - coordfn] [1] *sinv; 
y = coordfn] fo]*sinv + coordfn] fl]*cosv; 
z = coordfn] [2] ; 

fprintf (outfile, " %.81e, *.81e, *.81e\n",x,y,z); 

} 

} 


for(i=0; i<3; i++) { 

cosv = cos((double)i*PI*2./3. ) ; 
sinv = sin( (double) i*PI*2. /3. ) ; 

for(n=l; n<= N; n++) { 

x = coordfn] fO]*cosv - coordfn] [l]*sinv; 
y = coordfn] foj*sinv + coordfn] fl]*cosv; 
z = coordfn] f2] - h; 

fprintffoutfile, " £.81e, %.81e, *.81e\n",x,y,z) ; 

} 

} 

fclose(outfile) ; 

printf ("\nFile %s successfully written. \n" , outname) ; 


void initbeam(coord) 
double coordflOO] [3] ; 

{ • 

int i , n , j ; 

double x,y,z, xold, ksi,sqrt3, eps; 

sqrt3 = 3. ; 

sqrt3 = sqrt(sqrt3); 

x = b/(2.*sqrt3); 
y = b/2. ; 
coordfl] [0] = x; 
coordf 1] flj = y; 

coordfl] f2] = (x*x + y*y)/(4.*F); 

for(i=l; i<=M; i++) { 
n = i*(i + 1) + i; 
xold = x; 
eps = 1/2 . ; 

for( j=0; j<50; j++,eps/=2. ) { 
ksi = x + eps - xold; 
z = ksi*(ksi + 2.*xold)/(4.*F); 
if(ksi*ksi + z*z <- 1*1) 
x += eps; 

} 

coordfn] fO] = x; 
coordfn] fl] = y; 

coordfn] f2] = (x*x + y*y)/(4.*F); 

coordfn - 1][0] = (xold + x)/2.; 
coordfn - 1] fl] =0; 

coordfn - 1] [2] = x*x/(4.*F); . _ 
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} 


void fillsec(coord) 
double coord [100] [3] ; 

{ 

int i,n,p,q,r; 
double x,y,z; 

/* First fill in the other radial boundary 
*/ 

for(i=l; i<=M; i++) { 
n = i*(i + 1) + 1; 
sympoint (n, coord) ; 

} 

/* Then fill slanted rows and symmetrical points 
*/ 

for(i=l; i<=M; i++) { 
for(n=l; n<i; n++) { 

r = i*(i + 1) + 1 + 2*n; 

P = r - 2; 
q = r - 2*i; 

newpoint(p,q,r, coord) ; 
sympoint ( r , coord) ; 

} 

r = i*(i + 3) + 1; 
p = r-2; 
q = r - 1; 

newpo in t (p, q, r, coord) ; 

} 

} 


void sympoint (p, coord) 
int p; 

double coord[100] [3] ; 

{ 

double x,y,temp,cosv,sinv; 

temp = 2 . *PI/3 . ; 
cosv = cos (temp) ; 
sinv = sin(temp); 

x = coord[p] [0] ; 
y = coord[p][l]; 

coord[p + 1] [0] = x*cosv + y*sinv; 
coord [p + 1](1] = x*sinv - y#cosv; 
coordfp + 1] [2] = coord[p][2]; 


void newpoint(p*q,r, coord) 
int p,q,r; 

double coord [100] [3] ; 

{ 

int i,n; 
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double x,y,z,zold,Rbar[3] ,delR[3] , rhosq,d, temp.eps, mult; 


for(i=0,rhosq=0; i<3; i++) { 

Rbar[i] = (coord[p][i] + coord[q] [i] )/2; 
delR[i] = coord[p][i] - coord[q][i]; 
rhosq += delR[i]*delR[i] ; 


d = rhosq - delR[2]*delR[2] ; 

zold = (Rbar[0]*Rbar[0] + Rbar[ l]*Rbar [1] )/(4. *F) ; 
eps = zold/16; 

for(i=0,mult=l. ; i<50; i++,eps*=mult) { 
z = zold + eps - Rbar[2]; 
temp = 1*1 - rhosq*(0.25 + z*z/d) ; 
temp /= d; 
temp = sqrt(temp); 

x = Hbar[0] - z*delR[0]*delR[2]/d - delR[l]*temp; 
y = Rbar[l] - z*delR[lj*delR[2]/d + delR[0]*temp; 

if(x*x + y*y >= 4.*F*(zold + eps)) 
zold += eps; 
else mult = 0.5; 


} 

coordfr] [0] = x; 

coord[r][l] = y; 

coord[r] [2] = zold; 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

*/ 


GEN6STR.C - Generates the connectivity array for the six-spoke hybrid 
Pactruss. Input data is read from XYZ.DAT. 

The command is 

gen6str xyz 

If xyz is omitted, the default is SIXPAC. 

The connectivity array is created by the routines, is appended to 
the input nodal array and the combination is output to XYZ.DTA 
which is an ASCII file in the format needed for input to 
SEETRUSS.CMD. 

The types of struts are: 

Numbering 
External Internal 


1 

2 

3 

4 

5 

6 

7 

8 
9 


0 Non-folding Pactruss upper-surface struts 

1 Upper-surface Pactruss diagonals 

2 Upper-surface beam longerons and battens 

3 Upper-surface beam diagonals 

4 Non-folding Pactruss lower-surface struts 

5 Lower-surface Pactruss diagonals 

6 Lower-surface beam longerons and battens 

7 Lower-surface beam diagonals 

8 Pactruss verticals 

9 Beam verticals 

10 Core Pactruss ' diagonals 

11 Beam core diagonals 


John M. Hedgepeth 


Updated to MS C 5.0 


8/27/87 

12/29/87 


♦include <scitech.h> 
int M; 

int Ntypes = 12; 

int pacsurf ( ) , paccore ( ) , beamsurf ( ) , beamcore ( ) ; 


main(argc.argv) 
int argc; 
char. *argv [ ] ; 

{ 

int i , j , n , N, Nnodes , Nstruts , nl , n2 , type ; 

int Npacsurf , Npaccore , Nbeamsur f ^Nbeamcore ; 

int coni [100] [3] ,con2[100] [3] ,con3[100] [3] ,con4[100] [3] ; 

double x,y,z; 

char inname [ 50 ] , outname [ 50 ] ; 

FILE *infile,*outfile; 

strcpy ( inname , "SIXPAC . DAT" ) ; 
strcpy( outname, "SIXPAC. DTA") ; 

if(argc > 1) { 

strcpy ( inname , argv [ 1 ] ) ; 
strcpy(outname,argv[l] ) ; 
s treat (inname, ".DAT") ; 
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strcat( outname, " . DTA") ; 

} 

strupr( inname) ; 

strupr ( outname) ; 

if ( ( infile = fopen( inname, "rt")) == NULL) { 

printf ("\nCannot open %s. Aborted. \n", inname) ; 
exit(l) ; 

} 

if (fscanf( infile," *d, %lf, *lf, *lf, *lf *d", 

&M,&x,&x,&x,&x,8tNnodes) != 6) { 
printf("\nCannot read input parameters. Aborted. \n") ; 
exit(l) ; 


if((outfile = fopen(outname,"r”)) != NULL) { 
fflush(stdin) ; 

printf ("\nFile %s exists. Do you want to write over it? (Y/N) <N>”, 

outname) ; 

if(toupper(getch() ) != 'Y') 
exit(l) ; 

fclose(outfile) ; 

} 

outfile = fopen (outname, "wt") ; 

N = ( (M + 1)*(M + 2) )/2; 

Npacsurf = pacsurf(conl) ; 

Npaccore = paccore(con2) ; 

Nbeamsurf = beamsurf(con3) ; 

Nbeamcore = beamcore(con4) ; 

Nstruts = 6*(2*Npacsurf + Npaccore + 2*Nbeamsurf + Nbeamcore); 
fprintf (outfile, " Xd, Xd, ^d\n”, Nnodes, Nstruts, Ntypes ) ; 
for(i=0; i<Nnodes; i++) { 

if(fscanf( infile," *lf, *lf, *lf",&x,&y,&z) != 3) { 
printf ("Corrupt data for node *d. Aborted. \n" > i) ; 
exit(l); 

} 

fprintf (outfile," *.81e, %.81e, X.81e\n",x,y,z) ; 


fclose( infile) ; 

for(j=0; j<Npacsurf; j++) { 
for(i=0; i<6; i++) { 

nl = conl[j][0] + i*N; 
n2 = conl[j][l] + i*N; 

fprintf (outfile," Xd, %d, %d\n",nl,n2,conl[ j] [2] ) ; 
nl += 6*N; 
n2 + = 6W; 

fprintf (outfile," %d, %d, *d\n",nl,n2,conl[ j] [2] + 4); 

} 

} 

for(j=0; j<Npaccore; j++) { 
for(i=0; i<6; i++) { 

nl = con2[j][0J + i*N; 
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n2 = con2[ j ] [1] + i*N; 

fprintf (outfile, " *d, %d , %d\n",nl,n2,con2[ j] [2] ) ; 

} 

} 

for(j=0; j<Nbeamsurf ; j++) { 
for(i=0; i<6; i++) { 

nl = con3[j][0] + i*N; 
n2 = con3[j][l] + i*N; 
if(nl > 6*N) 
nl -= 6*N; 
if (n2 > 6*N) 
n2 -= 6*N; 

fprintf (outfile, " %d, %d, %d\n”,nl,n2,con3[ j] [2] ) ; 

nl += 6*N; 
n2 += 6*N; 

fprintf (outfile, " S»d, %d, %d\n”,nl,n2,con3[ j] [2] +4) 

} 

} 

for(j=0; j<Nbeamcore; j++) { 
for(i=0; i<6; i++) { 

nl = con4[j][0] + i*N; 
n2 B con4[j][l] + i*N; 
if (nl > 6*N) 
nl -= 6*N; 
if (n2 > 12*N) 
n2 -= 6*N; 

fprintf (outfile,” %d, %d, %d\n",nl,n2,con4[ j] [2] ) ; 

> 

} 

fclose(outfile) ; 

printf("\n\nData file *s successfully written.\n",outname) ; 


int pacsurf(buf) 
int buf [100] [3]; 

{ 

int i,n,p,nl,n2; 

for(n=l,p=0; n <= M; n++) { 
for(i=l; i<=n; i++) { 

nl = (n*(n + l))/2 + i + 1; 
n2 = ((n - l)*n)/2 + i; 

buf [p] [0] = nl; 
buf [p] [1] = nl - 1; 
buf [p] [2] = 1; 

p++; 

if (i == n) 
break; 

buf [p] [0] = nl; 
buf [pj [1] = n2; 
buf [p] [2] = 0; 

• 


A-14 



p++; 


buf [p] [0] = nl; 
buf [p] [1] = n2 + 1; 
buf [p] [2] = 0; 

p++; 

} 

} 

return p; 

} 


int paccore(buf) 
int buf [100] [3] ; 

{ 

int i,n,p,N,nl,n2; 

N = 6*( (M + 1)*(M + 2))/2; 

for(n=2,p=0; n<=M; n++) { 
for(i=0; i <= n; i++) { 

nl = (n*(n + l))/2 + i + 1; 


if(i>0 && i<n) { 
buf [p] [0] = nl; 
buf [p] [1] = nl + N; 
buf [p] [2] = 8; 

p++; 


if((n & 1) == 0) { 

n2 = ((n - l)*n)/2 + i; 

if(i>0 && i<n) { 
buf [p] [0] = nl; 
buf [p] [1] = n2 + N; 
buffpj [2] = 10; 

p++; 


buf [p] [0] = nl; 

buf [p] [1] = n2 + N + 1; 

buf [pj [2] = 10; 

p++; 


n2 = ((n + l)*(n + 2))/2 + i + 1; 

if(i > 0 && n < M) { 
buf [p] [0] = nl; 
buf [p] [1] = n2 + N; 
buf [p] [2] = 10; 

p++; 

} 

if(i < n && n < M) { 
buf [p] [0] = nl; 
buf[p][l] = n2 + N + 1; 
buf[p] [2] = 10; 

p++; 

} 
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} 


} 

return p; 


int beamsurf (buf ) 
int buf [100] [3] ; 

{ 

int i,nl,n2,p,N; 


N = 5*((M + 1)*(M + 2))/2; 


for(i=0,p=0; i<=M; i++) { 
nl = (i*(i + l))/2 + 1; 
n2 = ((i + l)*(i + 2) )/2 + 1; 

/* Beam battens 
*/ 

buf [p] [0] = nl; 
buf[p] [1] = nl + N + i; 
buf [p] [2] = 2; 

p++; 

if (i == M) 
break; 

/•* Beam longerons 
%/ 

buf [p] [0] = nl; 
buf [p] [1] = n2; 
buf [p] [2] = 2; 

p++; 

buf[p}[0] = nl + N + i; 
buf[p][l] = n2 + N + i + 1; 
buf [p] [2] = 2; 

p++; 

/% Beam diagonals 
*/ 

buf [p] [0] = nl; 

buf [p] [1] = n2 + N + i + 1; 

buf [p] [2] = 3; 

p++; 

} 

return p; 


int beamcore(buf) 
int buf [100] [3]; 

{ 


int i,nl,n2,p,N5,Nsurf; 

N5 = 5*( (M + 1)*(M + 2) )/2; 
Nsurf = 6*N5/5; 
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for(i=0,p=0; i<=M; i++) { 
nl = (i*(i + 1) )/2 + 1; 

buf [p] [0] = nl; 
buf[p][l] = nl + Nsurf; 
buf [p] [2] = 9; 

p++; 

if (i > 0) { 

buf [p] [0] = nl + N5 + i; 

buf [p] [ 1 ] = nl + N5 + i + Nsurf; 

buf [p] [2] = 9; 

p++; 

} 

buf [p] [0] = nl; 

buf[p][l] = nl + N5 + i + Nsurf; 
buf [p] [2] = 11; 

p++; 

if((i & 1) == 0) { 

n2n((i-l)*i)/2+l; 

if(i > 0) { 

buf [p] [0] = nl; 
buf[p][l] = n2 + Nsurf; 
buf [p] [2] = 11; 

p++; 

buf [p] [0] = nl + N5 + i; 
buf[p][l] = n2 + N5 + i - 1 + Nsurf; 
buf [p] [2] = 11; 

p++; 

} 

if (i < M) { 

n2 = ((i + l)*(i + 2) )/2 + 1; 

buf [p] [0] = nl; 

buf [p] [1] = n2 + Nsurf; 

buf [p] [2] = 11; 

p++; 

buf [p] [0] = nl + N5 + i; 
buf[pj[lj =n2+N5+i+l+ Nsurf; 
buf [p] [2] = 11; 

p++; 

} 

} 

} 

return p; 

} 
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/* GEN3STR.C - Generates the connectivity array for the three-spoke hybrid 
Pactruss. Input data is read from XYZ.DAT. The connectivity array 
is appended to the input nodal array and the combination is output 
to XYZ.DTA which is an ASCII file in the format needed for input to 
SEETRUSS.CMD. 


The types of struts are: 

Numbering 
External Internal 


1 

2 

3 

4 

5 

6 

7 

8 
9 


< 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Updated to MS C 5.0 


Non-folding Pactruss upper-surface struts 
Upper-surface Pactruss diagonals 
Upper-surface beam longerons and battens 
Upper-surface beam diagonals 
Non-folding Pactruss lower-surface struts 
Lower-surface Pactruss diagonals 
Lower-surface beam longerons and battens 
Lower-surface beam diagonals 
Pactruss verticals 
Beam verticals 
Core Pactruss diagonals 
Beam core diagonals 

John M. Hedgepeth 9/6/87 

12/29/87 


tinclude <scitech.h> 
int M; 

int Ntypes = 12; 


int pacsurf ( ) , paccore ( ) , beamsurf ( ) , beamcore ( ) ; 


main(argc,argv) 
int argc; 
char *argv[] ; 

{ 

int i, j,n,N,Nnodes,Nstruts,nl,n2,type; 

int Npacsurf , Npaccore , Nbeamsur f , Nbeamcore ; 

int coni [100] [3] ,con2[100] [3] ,con3[100] [3] ,con4[100] [3] ; 

double x,y,z; 

char inname [ 50 ] , out name [ 50 ] ; 

FILE *infile,*outfile; 

strcpy( inname, "TRIPAC.DAT") ; 
strcpy(outname, "TRIPAC.DTA") ; 


if(argc > 1) { 

strcpy( inname, argv[ 1] ) ; 
strcpy(outname,argv[l] ) ; 
s t rcat ( inname , " . D AT” ) ; 
strcat(outname,".DTA") ; 

} 


strupr ( inname) ; 
strupr(outname) ; 
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if( (infile = fopen(inname, "rt") ) == NULL) ( 

printf("\nCannot open Xs. Aborted. \n", inname) ; 
exit(l) ; 

} 

if(fscanf(infile, " Xd, %lf, %lf, 3>lf, %lf %d", 

&M,&x,&x,&x,&x,&Nnodes) != 6) { 
printf ("\nCannot read input parameters. Aborted. \n”) ; 
exit(l) ; 


if((outfile = fopen( out name, "r") ) != NULL) { 
fflush(stdin) ; 

printf ("\nFile 3»s exists. Do you want to write over it? (Y/N) <N>", 

out name) ; 

if(toupper(getch()) 1= *Y’) 
exit(l) ; 

fclose(outfile) ; 

} 

outfile = fopen(outname, "wt") ; 

N = M *(M + 3) +1; 

Npacsurf = pacsurf (coni) ; 

Npaccore = paccore(con2) ; 

Nbeamsurf = beamsurf (con3) ; 

Nbeamcore = beamcore(con4) ; 

Nstruts = 3*(2*Npacsurf + Npaccore + 2*Nbeamsurf + Nbeamcore); 
fprintf (outfile, " %d, %d , %d\n",Nnodes, Nstruts, Ntypes) ; 
for(i=0; i<Nnodes; i++) { 

if(fscanf( infile, " *lf, *lf, %lf",&x,&y,&z) != 3) { 
printf ("Corrupt data for node *d. Aborted. \n", i) ; 
exit(l); 

} 

fprintf (outfile," *.81e, *.81e, *.81e\n",x,y,z) ; 


fclose( infile) ; 

for(j=0; j<Npacsurf; j++) { 
for(i=0; i<3; i++) { 

nl = conl[j][0] + i*N; 
n2 = conl[j][lj + i*N; 

fprintf (outfile," Sid, *d, *d\n",nl,n2,conl[ j] [2] ) ; 
nl += 3*N; 
n2 += 3*N; 

fprintf (outfile, " *d, %d, \d\n",nl,n2,conl[ j] [2] +4); 

} 

} 

for(j=0; j<Npaccore; j++) { 
for(i=0; i<3; i++) { 

nl = con2[j](0] + i*N; 
n2 = con2(j][l] + i*N; 

fprintf (outfile," %d, Xd , %d\n",nl,n2,con2[ j][2] ) ; 

} 

} 


for(j=0; j<Nbeamsurf; j++) { 


A-19 


for(i=0; i <3 ; i++) { 

nl = con3( j] [0] + i*N; 
n2 = con3[j][l] + i*N; 
if (nl > 3*N) 
nl -= 3*N; 
if (n2 > 3*N) 
n2 -= 3*N; 

fprintf (outfile, " %d, %d, *d\n" , nl , n2 , con3 ( j] [2] ) ; 

nl += 3*N; 
n2 += 3*N; 

fprintf (outfile, " %d, %d, %d\n",nl,n2,con3[ j][2] + 4); 

} 

> 

for(j=0; j<Nbeamcore; j++) { 
for(i=0; i<3; i++) { 

nl = con4[j][0] + i*N; 
n2 = con4(j][l] + i*N; 
if (nl > 3*N) 
nl -= 3*N; 
if (n2 > 6*N) 
n2 -= 3*N; 

fprintf (outfile, ” %d, %d, *d\n",nl,n2,con4[ j] [2] ) ; 

} 

} 

fclose(outfile) ; 

printf ("\n\nData file %s successfully written. \n",outnaine) ; 


int pacsurf(buf) 
int buf [100] [3] ; 

{ 

int i,n,p,nl,n2; 

for(n=l,p=0; n <= M; n++) { 
for(i=l; i<=n; i++) { 

nl = n*(n + 1) + 2*i + 1; 
n2 = (n - l)*n + 2*i - 1; 

if(i < n) { 

buf [p] [0] = nl; 
buf [p] [1] = n2; 
buf [pj [2] = 1; 

p++; 

} 

buf [p] [0] = nl; 
buf [p] [1] = nl - 2; 
buf [p] [2] = 0; 


if (i == n) { 

buffp] [0] = nl; 
buf [p] [1] = nl - 1; 
buf [p] [2] = 0; 

p++; 
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buf [p] [0] = nl; 
buf[p][lj = nl - 2*(n + 1); 
buf[p] [2] = 1; 
p++; 

break; 

} 

buf [p] [0] = nl; 
buf [pj Cl] = n2 + 2; 
buf [p] [2] = 0; 

p++; 


buf [p] [0] = nl + 1; 
buf[ P ][l] = nl - 1; 
buf [p] [2] = 0; 

p++; 


buf [p] [0] = nl + 1; 
if(n == i + 1) 

buf [p] [1] = n2 + 2; 

else 

buf [p] [1] = n2 + 3; 
buf [p] [2] = 0; 

p++; 


buf [p] [0] = nl + 1; 
buf[p] [1] = n2 + 1; 
buf[p] [2] = 1; 

p++; 

} 

} 

return p; 

} 


int paccore(buf) 
int buf[100] [3] ; 

{ 

int i,n,p,N,nl,n2; * 

N = 3*(M*(M + 3) + 1); 

for(n=l,p=0; n<=M; n++) { 
for(i=0; i <= n; i++) { 

nl = n*(n + 1) + 2*i + 1; 

if(i>0) { 

buf [p] [0] = nl; 
buf [p] [1] = nl + N; 
buf [p] [2] = 8; 

p++; 


if(i != n) { 

buf [p] [0] = nl + 1; 
buf [p] [1] = nl + 1 + N; 
buf [p] [2] = 8; 

p++; 


} 
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if(((i + n) & 1) == 0) { 
if (i != n) { 

buf [p] [0] = nl; 

buf [p] [1] = nl + 2 + N; 

buf [p] [2] = 10; 

p++; 

buf [p] [0] = nl + 1; 
buf [p] [1] = nl + 3 + N; 
buf [p] [2] = 10; 

p++; 

if(i != 0) { 


buf [p] [0] 


nl; 



buf [pj [ 1] 


nl - 

2 + 

N; 

buffpj [2] 

= 

10; 



p++; 





buf [p] [0] 

= 

nl + 

l; 


buf [pj [1] 

= 

nl - 

l + 

N; 

buf [pj [2j 

= 

10; 



p++; 





buf [p] [0] 

= 

nl; 



buffpj [1] 


nl - 

2*n 

+ N; 

buf [pj [2 j 


10; 



p++; 





buf [p] [0] 

= 

nl + 

l; 


buf [p] [1] 


nl + 

l - 

2*n + N; 

buffpj [2 j 

= 

10; 



p++; 





if(n < M) 

{ 





buf[p] [0] = nl; 

buf [p] [1] = nl + 2*(n + 1) + N; 
buf [p] [2] = 10; 

p++; 

buf [p] [0] = nl + 1; 

buf[p][l] = nl + 1 + 2*(n + 1) + N; 

buf [p] [2] = 10; 

p++; 

} 

> 

} 

if (i == n) { 

buf [p] [0] = nl; 
buf[p][l] = nl - 1 + N; 
buf [p] [2] = 10; 

p++; 

buf[p] [0] = nl; 

buf [p] [1] = nl - 2 + N; 

buf [p] [2] = 10; 

p++; 


if(n < M) { 
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buf f p] [0] = nl; 
buf [p] [1] = nl 
buf [p] [2] = 10; 

p++; 


} 

} 

} 

return p; 


buf [p] [0] = nl; 
buf [p] [1] = nl 
buf [p] [2] = 10; 

p++; 


int beamsurf (buf) 
int buf [100] [3] ; 

{ 

int i,nl,n2,p,N; 

N = M*(M + 3) + 1; 


for(i=0,p=0; i<=M; i++) { 
nl = i*(i + 1) + 1; 
n2 = (i + l)*(i + 2) + 1; 

/* Beam battens 
*/ 

buf [p] [0] = nl; 
if ( i > 0) 

buf[p][l] = nl + 2*N + 1; 
else 

buf[p][l] = nl + 2*N; 
buf [p] [2] = 2; 

p++; 

if (i == M) 
break; 


/* Beam longerons 
*/ 

buf [p] [0] = nl; 
buf [p] [1] = n2; 
buf [p] [2] = 2; 

p++; 


if (i > 0) 

buf [p] [0] = nl + 2*N + 1; 

6XS6 

buf [p] [0] = nl + 2*N; 
buf[p][l] = n2 + 2*N + 1; 
buf [p] [2] = 2; 

p++; 


/* Beam surface diagonals 
*/ 


2*(n + 1) + 1 + N; 


2^(n + 1) + N; 


buf [p] [0] = nl; 
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buf [p] [ 1] = n2 + 2*N + 1; 
buffp] [2] = 3; 

p++; 


if (i > 0) 

buf[p] [0] = nl + 2*N + 1; 

else 

buf [p] [0] = nl + 2*N; 
buf [p] [1] = n2; 
buf [p] [2] =3; 

p++; 

} 

return p; 


int beamcore(buf ) 
int buf [ 100] [3] ; 

{ 

int i,nl,n2,p,N2, Nsurf; 

N2 = 2*(M*(M + 3) + 1); 
Nsurf = 3*N2/2; 

for(i=0,p=0; i<=M; i++) { 
nl = i*(i + 1) + 1; 


buffp] [0] = nl; 

buf [pj [1] = nl + Nsurf; . 

buf [pj [2] = 9; 

p++; 

if(i > 0) { 

buf [p] [0] = nl + N2 + 1; 

buf [p] [1] = nl + N2 + 1 + Nsurf; 

buf [p] [2] = 9; 

p++; 

} 

buf [p] [0] = nl; 
if(i > 0) 

buf [p] [1] = nl + N2 + 1 + Nsurf; 
gIsg 

buf [p] [1] = nl + N2 + Nsurf; 
buffp] [2] = 11; 

P++; 

if((i & 1) == 0) { 

n2 = (i - l)*i + 1; 

if(i > 0) { 

buf[p] [0] = nl; 

buf [pj [lj = n2 + Nsurf; 

buf [pj [2j = 11; 

p++; 

buf [p] [0] = nl + N2 + 1; 
buffp] [1] = n2 + N2 + 1 + Nsurf; 
buffp] [2] = 11; 
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i 


p++; 


} 

n2 = (i + l)*(i + 2) + 1; 

if(i < M) { 

buf [p] [0] = nl; 
buf[p][l] = n2 + Nsurf; 
buf [p] [2] = 11; 

p++; 

if (i > 0) 

buf[p][0] = nl + N2 + 1; 

else 

buf [p] [0] = nl + N2; 
buf [p] [ 1 J = n2 + N2 + 1 + Nsurf; 
buf [p] [2] = 11; 

p++; 

} 

} 

} 

return p; 
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TWORING. DTA 


73, 258, 12 

0 . 00000000e+000 , 0.00000000e+000, 4. 00000000e+000 
9 . 99999534e-001 , 5.77350000e-001, 8.33332556e-002 
2. 94193152e+000, 5.77350000e-001, 5.61768380e-001 
1 . 97096553e+000 , 2. 25911243e+000, 5.61768380e-001 
4 . 74468663e+000 , 5. 77350000e-001, 1. 4278365 le+000 
3 . 85313753e+000 , 2 . 22460999e+000 , 1 . 23722240e+000 
2 . 87234308e+000 , 3. 82034415e+000, 1.42783651e+000 
1. 12309793e-015, 1. 15470000e+000, 8.33332556e-002 
9. 70965992e-001, 2 . 83646243e+000 , 5.61768380e-001 
-9 . 70965992e-00 1 , 2 . 83646243e+000 , 5. 61768380e-001 
1 . 87234355e+000 , 4. 397694 15e+000, 1 . 42783651e+000 
6. 11750234e-015, 4.44921998e+000, 1.23722240e+000 
-1 . 87234355e+000 , 4. 39769415e+000, 1.4278365 le+000 
-9 . 99999534e-00 1 , 5.77350000e-001, 8. 33332556e-002 
-1 . 97096553e+000 , 2.25911243e+000, 5.61768380e-001 
-2.94193 152e+000 , 5.77350000e-001, 5. 61768380e-001 
-2 . 87234308e+000 , 3.82034415e+000, 1.4278365 le+000 
-3 . 85313753e+000 , 2 . 22460999e+000 , 1.23722240e+000 
-4 . 74468663e+000 , 5.77350000e-001, 1.4278365 le+000 
-9 . 99999534e-00 1 , -5.77350000e-001, 8.33332556e-002 
-2.94193152 e+000, -5. 77350000e-001, 5. 61768380e-001 
-1 . 97096553e+000 , -2. 25911243e+000, 5 . 61768380e-001 
-4 . 74468663e+000 , -5.77350000e-001, 1.4278365 le+000 
-3 . 853 13753e+000 , -2 . 22460999e+000 , . 1 . 23722240e+000 
-2 . 87234308e+000 , -3.82034415e+000, 1.4278365 le+000 
-4 . 57397792e-015 , -1 . 15470000e+000 , 8 . 33332556e-002 
-9. 70965992e-001 , -2.83646243e+000, 5 . 61768380e-001 
9. 70965992e-001 , -2.83646243e+000, 5. 61768380e-001 
-1 . 87234355e+000 , -4. 397694 15e+000, 1.42783651e+000 
-1.94142662e-014, -4.44921998e+000, 1 . 23722240e+000 
1 . 87234355e+000 , -4. 39769415e+000, 1.4278365 le+000 
9 . 99999534e-00 1 , -5.77350000e-001, 8.33332556e-002 
1.97096553e+000, -2.25911243e+000, 5.61768380e-001 
2.94193152e+000, -5.77350000e-001, 5.61768380e-001 
2 . 87234308e+000 , -3. 820344 15e+000, 1.4278365 le+000 
3.85313753e+000, -2 . 22460999e+000 , 1.23722240e+000 
4 . 74468663e+000 , -5.77350000e-001, 1.4278365 le+000 
9 . 99999534e-00 1 , 5.77350000e-001, -1.91666674e+000 
2.94193152e+000, 5. 77350000e-001 , -1.43823162e+000 
1 . 97096553e+000 , 2.25911243e+000, -1. 43823 162e+000 
4 . 74468663e+000 , 5.77350000e-001, -5.72163486e-001 
3. 85313753e+000, 2 . 22460999e+000 , -7.62777599e-001 
2 . 87234308e+000 , 3 . 82034415e+000 , -5.72163486e-001 
1. 12309793e-015, 1.15470 000e+000, -1. 91666674e+000 
9. 70965992e-001 , 2 . 83646243e+000 , -1.43823162e+000 
-9 . 70965992e-00 1 , 2.83646243e+000, -1.43823162e+000 
1 . 87234355e+000 , 4. 397694 15e+000, -5.72163486e-001 
6. 11750234e-015, 4.44921998e+000, -7.62777599e-001 
-1 . 87234355e+000 , 4.39769415e+000, -5. 72163486e-001 
-9. 99999534e-001 , 5.77350000e-001, -1.91666674e+000 
- 1 . 97096553e+000 , 2. 25911243e+000, -1. 43823 162e+000 
-2 . 94193152e+000 , 5.77350000e-001, -1.43823162e+000 
-2 . 87234308e+000 , 3. 820344 15e+000, -5.72163486e-001 
-3 . 85313753e+000 , 2 . 22460999e+000 , -7.62777599e-001 
-4 . 74468663e+000 , 5.77350000e-001, -5.72163486e-001 
-9 . 99999534e-00 1 , -5.77350000e-001, -1.91666674e+000 
-2.94193152e+000, -5.77350000e-001, -1.43823162e+000 
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-1 . 97096553e+000 , -2. 25911243e+000, -1.43823162e+000 
-4 . 74468663e+000 , -5.77350000e-001, -5.72163486e-001 
-3. 85313753e+000, -2. 22460999e+000 , -7.62777599e-001 
-2 . 87234308e+000 , -3.82034415e+000, -5. 72163486e-001 
-4 . 57397792e-0 15 , -1. 15470000e+000, -1.91666674e+000 
-9.70965992e-001, -2.83646243e+000, -1.43823162e+000 
9. 70965992e-001 , -2 . 83646243e+000 , -1.43823162e+000 
-1 . 87234355e+000 , -4.39769415e+000, -5. 72163486e-001 
-1 . 94142662e-014 , -4.44921998e+000, -7. 62777599e-001 
1 . 87234355e+000 , -4.39769415e+000, -5.72163486e-001 
9. 99999534e-001 , -5.77350000e-001, -1.91666674e+000 

I . 97096553e+000 , -2.25911243e+000, -1.43823162e+000 
2.94193152e+000, -5.77350000e-001, -1.43823162e+000 
2 . 87234308e+000 , -3. 820344 15e+000, -5.72163486e-001 
3 . 853 13753e+000 , -2. 22460999e+000, -7.62777599e-001 
4 . 74468663e+000 , -5.77350000e-001, -5.72163486e-001 
3, 2, 1 
39, 38, 5 
9, 8, 1 
45, 44, 5 
15, 14, 1 
51, 50, 5 
21 , 20 , 1 
57, 56, 5 
27, 26, 1 
63, 62, 5 
33, 32, 1 
69, 68, 5 
5, 4, 1 
41, 40, 5 

II, 10, 1 
47, 46, 5 
17, 16, 1 
53, 52, 5 
23, 22, 1 
59, 58, 5 
29, 28, 1 
65, 64, 5 
35, 34, 1 
71, 70, 5 
5, 2, 0 
41, 38, 4 
11, 8, 0 
47, 44, 4 
17, 14, 0 
53, 50, 4 
23, 20, 0 
59, 56, 4 
29, 26, 0 
65, 62, 4 
35, 32, 0 
71, 68, 4 
5, 3, 0 
41, 39, 4 
11, 9, 0 
47, 45, 4 
17, 15, 0 
53, 51, 4 
23, 21, 0 
59, 57, 4 
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1, 31, 2 
37, 67, 6 
7, 1, 2 
43, 37, 6 
13, 7, 2 
49, 43, 6 
19, 13, 2 
55, 49, 6 
25, 19, 2 
61, 55, 6 
31, 25, 2 
67, 61, 6 

I 1, 2, 2 

37, 38, 6 
7, 8, 2 
43, 44, 6 
13, 14 ", 2 
49, 50, 6 
19, 20, 2 

i 55, 56, 6 

25, 26, 2 
61, 62, 6 
31, 32, 2 
67, 68, 6 
31, 33, 2 
67, 69, 6 
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1, 3, 2 
37, 39, 6 
7, 9, 2 
43, 45, 6 
13, 15, 2 
49, 51, 6 
19, 21, 2 
55, 57, 6 
25, 27, 2 
61, 63, 6 

1, 33, 3 

37, 69, 7 

7, 3, 3 

43, 39, 7 

13, 9, 3 

49, 45, 7 

19, 15, 3 

55, 51, 7 

25, 21, 3 

61, 57, 7 

31, 27, 3 

67, 63, 7 

2, 33, 2 

38, 69, 6 

8, 3, 2 

44, 39, 6 

14, 9, 2 

50, 45, 6 

20, 15, 2 

56, 51, 6 

26, 21, 2 

62, 57, 6 

32, 27, 2 

68, 63, 6 

2, 4, 2 

38, 40, 6 

8 , 10 , 2 

44, 46, 6 

14, 16, 2 

50, 52, 6 

20 , 22 , 2 

56, 58, 6 
26, 28, 2 

62, 64, 6 

32, 34, 2 

68, 70, 6 

33, 36, 2 

69, 72, 6 

3, 6, 2 

39, 42, 6 

9, 12, 2 

45, 48, 6 

15, 18, 2 

51, 54, 6 

21, 24, 2 

57, 60, 6 
27, 30, 2 

63, 66, 6 
2, 36, 3 
38, 72, 7 
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/* SEE2PAL.C - Reads geometry from XYZ.DTA (SEETRUSS format), strut cross 

* section and material properties for each type of strut from file 

* XYZ.PRP. Produces file XYZ.TXT which is suitable for input into 

* the MSC PAL2 finite element analysis program. 

* 

* The file XYZ.PRP has a row, for each type, of the following quanti- 

* ties separated by white space for each type: type number (starting 

* at 1), strut area A, Young’s modulus E, material density rho, and 

* the strut Euler buckling stress S_euler. The material density is 

* adjusted to include joint and payload masses. 

* 

* John M. Hedgepeth 8/30/87 

* 

* Updated to MS C 5.0 12/30/87 

* Added thermal properties to material specs 5/11/88 

* 

* 

*/ 

#include <scitech.h> 

char *getdataline(FILE *); 

main (argc, argv) 
int argc; 
char *argv [ ] ; 

{ 

int i , n , N, Ntypes , Nnodes , Nstruts , nl , n2 , type , *p_strut , *ptr ; 
double x,y,z, A,E, rho,S_euler; 

char dtaname[30] ,prpname[30] , outname[30] ,name[30] , *line; 

FILE *dtafile,*prpfile,*outfile; 

strcpy(dtaname, "SIXPAC.DTA”) ; 
s t rcpy ( prpname , "SIXPAC.PRP”) ; 
strcpy(outname, "SIXPAC.TXT”) ; 
strcpy(name, "SIXPAC”) ; 

if(argc > 1) { 

s t rcpy ( dtaname , argv [ 1 ] ) ; 
strcpy(prpname,argv[l] ) ; 
s t rcpy ( outname , argv [ 1 ]) ; 
s trcpy( name, argv [1] ) ; 

strcat(dtaname,”.DTA”) ; 
s treat (prpname, ” . PRP" ) ; 
streat ( outname TXT” ) ; 

} 

strupr(dtaname) ; 
strupr(prpname) ; 
strupr( outname) ; 
strupr(name) ; 

if((dtafile = fopen(dtaname, "rt”) ) == NULL) { 

printf (”\nCannot open %s. Aborted. \n”, dtaname) ; 
exit(l); 

} 

if ( (prpfile = fopen( prpname, "rt”) ) == NULL) { 

printf(”\nCannot open *s. Abor ted. \n", prpname ) ; 
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} 


exit( 1) ; 


line = getdataline(dtafile) ; 

if(sscanf(line," %d, %d, *d",&Nnodes,&Nstruts,&Ntypes) != 3) { 

printf("\nCannot read input parameters from %s. Aborted. \n”, dtaname) ; 
exit(l) ; 

} 

/* Allocate space for strut connection array 
*/ 

if((p_strut = (int *)calloc(3*Nstruts,sizeof (int) ) ) == NULL) { 

printf ("\nCannot allocate space for strut array. Aborted. \n") ; 
exit(l) ; 

} 

if((outfile = fopen( out name, "r”) ) ! = NULL) { 
fflush(stdin) ; 

printf ("\nFile %s exists. Do you want to write over it? (Y/N) <N>” , 

outname) ; 

if(toupper(getch() ) != *Y*) 
exit(l); 

fclose(outfile) ; 

} 

if( (outfile = fopen( outname, "wt") ) == NULL) { 

printf ("\Cannot open %s. Aborted. \n", outname) ; 
exit(l); 


/% Start with truss name 
*/ 

fprintf (outfile, "TITLE %s TRUSS\n\n'\name) ; 

/* Discard first (zeroeth) node 

*/ 

line = getdataline(dtafile) ; 

if(sscanf(line,'* *lf, ^lf, %lf",&x,&y,&z) != 3) { 
printf("\nCannot scan node 0. Aborted. \n") ; 
exit(l) ; 

} 


/* Transfer node coordinates 
*/ 

fprintf (outfile, "NODAL POINT LOCATIONS l\n"); 

for(n=l; n<Nnodes; n++) { 

line = getdataline(dtafile) ; 

if(sscanf(line,” *lf, S»lf, \lf ,, ,Stx,&y,&z) != 3) { 
printf ("\nCannot scan node %d. Aborted. \n",n) ; 
exit(l) ; 

} 

fprintf (outfile, "*-3d *12.61f, %12.61f, %12.61f\n",n,x,y,z) ; 

} 

fprintf (outfile, "\n") ; 

/* Discard rotation degrees of freedom 
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*/ 

fprintf (outfile, "ZERO l\nRA OF ALL\n\n") ; 

/* Select type of mass lumping 
*/ 

fprintf (outfile, "LUMPED MASS TYPE l\n"); 

/* Fill the strut connection array 
*/ 

ptr = p_strut; 

for(i=0; i<Nstruts; i++) { 

line = getdataline(dtafile) ; 

if(sscanf(line," 3»d, %d, %d",&nl,&n2,&type) != 3) { 

printf ("\nCannot read data for strut %d. Abort. \n", i) ; 
exit(l) ; 

} 

. *ptr++ = type; 

*ptr++ = nl; 

*ptr++ = n2; 

} 

fclose(dtafile) ; 

/* Transfer properties, connectivity, and ANALYZE for each type 
*/ 

fprintf (outfile, "C\nC ********** Struts ****************\nC\n'') ; 
fprintf (outfile, "\nANALYZE 0\n") ; 

N = Nnodes/12; 

for(n=0; n<Ntypes; n++) { 

line = getdataline(prpfile) ; 

if (sscanf(line, ” *d %lf %lf *lf %lf",&type,&A,&E,&rho,&S_euler) 

!= 5) { 

print f("\nCannot read properties for type %d. Aborted. ", n+1) ; 
exit(l) 

} 

} 

fprintf (outfile, M \nC ** Type *d **\n\n",type) ; 

fprintf (outfile, "MATERIAL PROPERTIES *1E, 0.0, *lf, 0.3, %1E, * 

"1.0, 0.0\n", E,rho,S_euler) ; 
fprintf (outfile, "BEAM TYPE 1, *lE\n",A) ; 

ptr = p_strut; 
for(i=0; KNstruts; i++) { 
if(*ptr++ == n) { 
nl = *ptr++; 
n2 = *ptr++; 

fprintf (outfile, "CONNECT %d TO *d\n",nl,n2) ; 

} 

else { 

ptr++; 

ptr++; 

} 

} 

} 

fprintf (outfile, "\END DEFINITION\n") ; 
fclose(prpfile) ; . 



fclose(outfile) ; 

printf ("\nPal2 input file %s successfully written. \n",outname) 

} 

/ '* ******************************************************** 

*/ 

char *getdataline(FILE *infile) 

{ 

int chr,old; 
char *ptr; 

static char line[81]; 

while(fgets(line, 80, infile) != NULL) { 
ptr = line; 

while(isspace( (chr = (int)*ptr))) 
ptr++; 

if(isdigit(chr) !! chr == !! chr == *+’ ! ! chr == ’-’) 

return ptr; 

} 

if (feof(infile) ) { 

printf ("\007\nEnd of file reached. \n") ; 
return line; 

} 

else { 

printf (”\007\nError in reading pal2 file. Abort. \n"); 
exit(l) ; 

} 

} 

/******************************************************** 

*/ 
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/* BESTFIT.C - Determines the location x0,y0,z0 of the vertex and the 

* orientation angles phix.phiy of the best-fit paraboloid for 

* a deformed surface. The focal length of the paraboloid is held 

* constant. The parameters of the reference paraboloid and the 

* coordinates of nodes on the undeformed surface are read from the 

* file XYZ.DAT which was formed by the command 

* 

* genNpac.exe xyz (N = 6 or 3) 

* 

* The values of the deformations at each point are read from the 

* file ABC.DEF which is obtained from the output of MSC Pal2 STAT2 

* by using a text editor to extract the title and the list of dis- 

* placements . 

* 

* The method used is to find the values of xO,yO,zO,phix,phiy which 

* minimize the sum of (Un cosalfa)*(Un cosalfa)over all the nodes, 

* where 

* 

* Un cosalfa = {uz - zO + x*phiy -y*phix 

* - (ux - xO - z*phiy)*x/(2*F) 

* - (uy - yO + z*phix)*y/(2*F)}/{l + z/(2*F)} 

* 

* is one-half of the error in the reflected path length for the 

* displaced paraboloid. A weighting factor of cosalfa is used. 

* The displacements xO and yO are constrained to be related to the 

* rotations so that the focal point is stationary. Thus 

* 

* xO = -F*phiy 

* yO = F*phix 

* 

* Results are reported to the screen as values of the paraboloid 

* displacement and rotation as well as the rms value of the residual 

* half-path-length error. 

* 

* j.m.hedgepeth 9/14/87 

* 

* Updated for MS C 5.0 12/30/87 

.*/ 

finclude <scitech.h> 
int M; 

double disp[5],F; 

main() 

{ 

int i, j,m,flag, lpa[5] ; 

double x,y,z,*p_crd,*p_def,*p_real,c[6] ,d[6] ,Mat[3] [4] ,w,W; 

FILE *datfile,*deffile; 

char datname[50] ,defname[50] , linebuf [100] , *p_chr,*n; 
char title[80] ; 

printf("\n\nBESTFIT reads input from xyz. DAT and abc.DEF.\n") ; 
printf ("Enter pathname of .DAT file? "); 
gets(datname) ; 
putchar( V\n' ) ; 

if((n = strchr(datname, ’ . * ) ) == NULL) 
s treat ( datname , " . DAT" ) ; 
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if((datfile = fopen(datname, "rt") ) == NULL) { 

printf("\nCannot open %s for reading. Aborted. \n’\ datname) ; 
exit(l) ; 

} 

printf ("Enter pathname of .DEF file? "); 
gets(defname) ; 
putchar( ’ \n’ ) ; 

if((n = strchr(defname, ’ . ’ )) == NULL) 
strcat(defname, ".DEF") ; 

if ( (deffile = f open (def name, "rt") ) == NULL) { 

printf("\nCannot open *s for reading. Aborted. \n",defname) ; 
exit(l) ; 

} 

if (fscanf (datfile, " *d, *lf, %lf, %lf, *lf",&i,&F,&x,&y,&z) != 5) { 
printf ("\nCannot read parameters from .DAT file. Abort. \n"); 
exit(l) ; 

} 

if (fscanf (datfile," *d",&M) != 1) { 

printf("\nCannot read number of nodal points. Abort. \n"); 
exit(l) ; 

} 

/* Ignore first node (focal point) and divide M - 1 by 2 to use only upper 
* surface points. 

*/ 

M = (M - l)/2; 

if (fscanf (datfile," *lf, *lf, *lf",&x,&y,&z) != 3) { 

printf ("\nCannot scan coordinates of node 0. Abort. \n"); 
exit(l) ; 

} 

/* Allocate space for arrays 
*/ 

if((p_crd = (double *)calloc(3*M,sizeof (double))) == NULL) { 

printf("\nCannot allocate memory for coordinate array. Abort. \n"); 
exit(l) ; 

} 

if ( (P_def = (double *)calloc(3*M,sizeof (double))) == NULL) { 

printf(”\nCannot allocate memory for deflection array. Abort. \n"); 
exit(l) ; 

> 

/* Fill coordinate array 
*/ 

p_real = p_crd; 
for(i=0; i<M; i++) { 

if (fscanf (datfile," *lf, *lf, %lf",&x,&y,&z) != 3) { 

printf ("\nCannot scan coordinates for node %d. Abort. \n", i+1) ; 
exit(l) ; 

} 

*p_real++ = x; 

*p_real++ = y; 

*p_real++ = z; 

} 

fclose(datfile) ; 
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/* Examine first lines of .DEF file until we reach deflection data. Pick up 
* title on the way. 

*/ 

/* Search each line for the title or for start of deflections 
*/ 


flag = 1; 
while(flag) { 

fgets(linebuf , 100, deffile) ; 

if (sscanf ( linebuf , " Xd *lf *lf *lf",&j,&x,&y,&z) == 4) 
break; 

p_chr = linebuf; 

while(*p_chr == ’ ’ !! *p_chr == ’\t’ I! *p_chr == *\n’) 
p_chr++; 
if(*P_chr) { 

strcpy( title, p_chr) ; 
flag = 0; 

} 

} 

if (flag) 

strcpy(title," TRUSS"); 

@Xs6 { 

while(fgets(linebuf, 100, deffile) != NULL) { 

if (sscanf (linebuf," Xd Xlf Xlf *lf”,&j,&x,&y,&z) == 4) { 
flag = 1; 
break; 

} 

} 

} 

if (flag == 0) { 

printf(”\nError reading start of .DEF file. Abort. \n"); 
exit(l) ; 


/* Fill deflection array 
*/ 


p_real = p_def; 


for(i=0; i< j— 1 ; i++) { 
*p_real++ = 0; 
*p_real++ = 0; 
*p_real++ = 0; 


*p_real++ = x; 
*p_real++ = y; 
*p_real++ = z; 

i++; 


for(; i<M; i++) { 

if (fgets( linebuf, 100, deffile) == NULL) { 

printf ("\nCannot read .DEF file for node *d. Abort. \n",i + 1); 
exit(l) ; 

} 

if (sscanf (linebuf," Xd Xlf *lf *lf",&j,&x,&y,&z) != 4) { 

printf("\nCannot scan deflections of node %d. Abort. \n",i + 1); 
exit(l) ; fl ,o 



} 

for( ; i< j-1; i++) { 
*p_real++ = 0; 
*p_real++ = 0; 
*p_real++ = 0; 


*p_real++ = x; 
*p_real++ = y; 
*p_real++ = z; 


fclose(deffile) ; 

/* Prepare to fill matrix 
*/ 


for(i=0; i<4; i++) { 
for( j=0; j<3; j++) 

Mat[ j] [i] = 0; 
d[i] = 0; 

} 

W = 0; 

for(m=0; m<M; m++) { 
x = *p_crd++; 
y = *p_crd++; 
z = *p_crd++; 

x /= 2.*F; 
y /= 2.*F; 
z /= 2.*F; 

w = l./(l. + 2.*z); 

c[0] = x*w; 
cfrj = y*w; 

c[2] = -w; 
c[3] = -y; 
c[4] = x; 

w = sqrt(w) ; 

x = *p_def++; 
y = *p_def++; 
z = *p_def++; 

c[5] = c[0]*x +c[l]*y + c[2]*z; 

for(i=0; i<4; i++) { 
for( j=0; j<3; J++) 

Mat[j][i] += c[i + 2]*c[j + 2]*w; 
d[i] += c[i + 2]*c[5]*w; 

} 

W += w; 


for(i=0; i <4 ; i++) 
d[i] /= W; • 
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w = .000001; 

dcrout(4,3,l,Mat,w,&W,&i, lpa) ; 
W = d[3] ; 

for(i=0; i<3; i++) { 

disp[i] = Mat[i] [3] ; 

W -= d[i]*disp[i] ; 


for(i=0; i<25; i++) 
putchar( ’ \n* ) ; 

printf (" %s\n\n", title) ; 

printf("\tRms distortion from best-fit paraboloid 
printf ("\tVertical displacement of paraboloid: zO 

printf("\tRotation of paraboloid about: x axis 

printf ("\t y axis 

for(i=0; i< 10 ; i++) 
putchar( ’\n’ ) ; 

fflush(stdin) ; 
getch( ) ; 


%le\n\n" , sqrt ( W) ) ; 
*le\n\n" , disp [ 0 ] ) ; 
%le\n”,disp[l]/2./F) ; 
%le\n",disp[2]/2./F) ; 
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